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ABSTRACT 


Total  time  on  test  (TTT)  plots  provide  a useful  graphical 
method  for  tentative  identification  of  failure  distribution 
models.  Identification  is  based  on  properties  of  the  TTT 
transform.  New  properties  of  the  TTT  transform  disti  Lbution 
are  obtained.  In  particular,  it  is  shown  that  a non-IFRA 
distribution  may  have  an  anti-starshaped  transform.  Hence, 
TTT  transforms  may  only  be  useful  for  determining  local 
properties  of  the  failure  rate  function  and  not  the  failure 
rate  average  function. 
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GEOMETRY  OF  THE  TOTAL  TIME  ON  TEST  TRANSFORM 

by 

Richard  E . Barlow 

1.  INTRODUCTION 

The  geometry  of  the  total  time  on  test  transform  is  helpful  in  inter- 
preting total  time  on  test  data  plots  [cf.  Barlow  and  Campo  (1975)].  In 
particular,  it  is  possible  to  infer  tentative  probability  distribution 
models  based  on  total  time  on  test  plots. 

Let  F be  a failure  distribution,  i.e.,  F(0  ) = 0 and  F = 1 - F . 

Define 

F~l(t> 

H^1(t)  = J F(x)dx  0 <_  t <_  1 , 

0 


the  total  time  on  test  transform  of  F . It  is  easy  to  verify  that,  H^,  , 
the  inverse  of  is  a distribution  function.  Also,  H has  support 

in  [0,0]  if  0 is  the  mean  of  F , since 


F"i(l) 


J 


F(x)dx  = J xdF(x)  = 0 
0 


by  an  integration  by  parts.  It  is  easy  to  verify  that  if  F(x)  = 1 - e 
then  the  corresponding  HF(x)  = x/0  for  0 _<  x £ 0 . The  result  that 
our  transform  carries  the  exponential  distribution  into  the  rectangular 
distribution  on  [0,0]  is  important.  It  is  easy  to  verify  that  in  this 
case  Hj,  < F ; i.e.,  F ^Hp(x)  is  convex  for  0 <_  x <_  6 . 
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If  F(x)  = x for  0 <_  x 1 ; i.e.,  F is  uniform  on  [0,1]  , 

then  Hj,(x)  = 1 - /I  - 2x  for  0 j<  x <_  h , since  F has  mean  % 
in  this  case.  Again  < F . In  general  Hp^F(x)  * J* F(u)du  is 

° _i 

concave  since  F is  nonincreasing.  Hence,  the  inverse  F Hy(x)  is 
convex  for  0 <_  x <_  9 and 

HF  i F 

for  all  F with  F(0  ) = 0 . 

As  was  proved  in  Barlow,  Bartholomew,  Bremner  and  Brunk  (1972), 
total  time  on  test  data  plots  tend  to  the  total  time  on  test  transform 
of  the  underlying  failure  distribution  as  the  sample  size  tends  to  infinity. 
In  order  to  interpret  total  time  on  test  data  plots,  we  need  to  understand 
the  relationship  between  F and  its  transform.  The  following  table 
summarizes  the  connections. 


Life  Distribution 

Hazard  Function 

Total  Time  on  Test 
Transform  Distribution 

F 

R = log  F 

«F 

Exponential 

linear 

linear 

IFR 

convex 

<*=£> 

convex 

DFR 

— 

concave 

concave 

IFRA 

starshaped 

- 

starshaped 

DFRA 

<*<> 

anti-star shaped 

■fr 

anti-star shaped 

TABLE  1 

Logical  Connections  Between  Life  Distributions, 

Hazard  Functions  and  TTT  Transform  Distributions 

__ ^ 


3 


o (x) 

A function  g defined  on  [0,b)  such  that  ' is  nondecreasing 
on  [0,b)  is  said  to  be  starshaped  with  respect  to  the  origin.  If 
G(x)  = 1 - e , then  F is  IFRA  (for  increasing  failure  rate  average 

g-^Ffx)  —1 

if  and  only  if  is  nondecreasing  for  0 <_  x <_  F (1)  . The 

function  G "^(x)  = R(x)  is  said  to  be  starshaped  with  respect  to  the 

origin.  As  the  last  two  implications  indicate,  IFRA  and  DFRA  distribution 

families  are  not  characterized  by  corresponding  properties  of  the  TTT 

transform  distribution.  However  IFR  and  DFR  distribution  families  are 

characterized  by  corresponding  properties  of  the  TTT  transform  distribution. 

To  verify  the  implications  in  the  table  for  the  IFR  (DFR)  case, 
first  assume  F absolutely  continuous  with  failure  rate  function,  r . 

If  F is  IFR  (DFR),  then 

is  decreasing  (increasing)  in  x which  implies  is  concave  (convex) ; 

i.e.,  Hp,  is  convex  (concave).  Conversely,  if  Hp/  is  concave  (convex), 
the  failure  rate  function  is  increasing  (decreasing).  To  see  this,  note 
that  every  IFR  (DFR)  distribution  can  be  approximated  arbitrarily  closely 
by  an  absolutely  continuous  IFR  (DFR)  distribution.  Since  the  limit  of 
a sequence  of  concave  (convex)  transforms  is  concave  (convex)  on  [0,1]  , 
it  follows  that  F is  IFR  (DFR)  if  and  only  if  Hp,^  is  concave  (convex) . 

The  IFRA  distributions  govern  the  lifelength  of  coherent  systems 
with  statistically  independent  components  whose  life  distributions  are 
IFR  (or,  more  generally,  IFRA).  [Birnbaum,  Esary  and  Marshall  (1966)  or 
Barlow  and  Proschan  (1975)].  They  also  arise  in  other  reliability  contexts. 


= 1 

t=F (x)  ' r(x) 


For  these  reasons,  we  are  interested  in  the  transforms  of  IFRA  distributions 


In  the  next  section  we  show  that  if  F is  IFRA,  then  its  transform  dis- 

Vx> 

tribution,  , is  starshaped;  i.e.,  — - — is  nondecreasing  in 

0 < x < 0 . Unfortunately,  the  converse  is  not  true. 
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2.  PRESERVATION  OF  PARTIAL  ORDERINGS  ON  CLASSES  OF  FAILURE  DISTRIBUTIONS 

Let  R(x)  = -log  F(x)  be  the  hazard  function  of  F as  before  and 
let  G(x)  = 1 - e x . Observe  that  G ^F(x)  = R(x)  so  that  if  F is 
IFR,  G ^F(x)  is  convex  on  the  support  of  F and  conversely.  If  F 
G~^f (x) 

is  IFRA,  is  nondecreaslng  in  x ^ 0 and  conversely.  This 

leads  to  a partial  ordering  on  the  space  of  failure  distributions  which 
we  call  "star  ordering."  Let  F be  the  class  of  continuous  distributions 
on  [0,®)  and  {deg.}  , the  class  of  degenerate  distributions. 


Definition: 

F1  ^ F2  (1,e’>  *s  star  ordered  with  respect  to  F£  if  F^  , 

F2lpi(x) 

F2  e FU(deg.}  and  is  nondecreasing  in  x for  0 

According  to  this  definition,  every  distribution  in  F is  star 

_x« 

ordered  with  respect  to  a degenerate  distribution.  Let  F^Cx)  = 1 - e 

for  x ^ 0 and  a > 0 . It  is  easy  to  show  that  if  0 < a^  < 02  , then 

F < F . Since  F has  failure  rate  ax  , it  is  clear  that  the 

02  * a^  a * 

failure  rate  of  F is  "increasing  faster"  than  the  failure  rate  of  F 

a2  al 

If  0 < a < 1 , F is  DFR.  If  a > 1 , F is  IFR  and  F,  is  exponential. 

a ’a  1 r 

Definition: 

F^  < F2  (i.e.,  F^  is  convex  ordered  with  respect  to  F2  if 

-1  -1  \ 

Fi  , Fj  e F U {deg.}  and  F2  F^(x)  is  convex  in  x for  0 x F^  (1) J . 

It  is  not  hard  to  show  that  c-ordering  implies  star  ordering,  but 

not  conversely.  Our  main  theorem  is  that  the  TTT  transform  distribution 

preserves  both  orderings. 


< x < 


F^d) 
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Theorem  2.1: 

Let  , F2  e F . 

(a)  If  Fx  < F2  , than  . 

«.)  If  Fx  < F2  , than  iIFi  J “f2  • 

The  following  corollary  provides  the  primary  application  of  the 
theorem. 

Corollary  2.2: 

If  Fj^  , F2  £ F and  F^  < F2  , then 

»;V> 

(a)  — = is  nondecreasing  in  0 < t < 1 ; 

tfw  - ~ 

H^(t)  H^(t) 

(b)  — y >_  — y for  0 <_  t <_  1 . 

Hi  (1)  C(l) 

F1  T2 

Proof  of  Corollary; 

h;\  (x> 

By  Theorem  2.1,  Part  (b),  H < H so  — - — is 

F1  F2  x 

nondecreasing  in  0 <_  x _<  F^d)  . Let  t * Hy  (x)  and  Part  (a)  of 
the  corollary  is  immediate,  (b)  is  a trivial  consequence  of  (a) . | | 

Figures  1 and  2 are  graphical  plots  of  the  scaled  transforms  of 


gamma  and  Veibull  distributions.  They  visually  confirm  Part  (b)  of 
the  corollary. 


10 


Proof  of  Theorem  2.1: 


(a)  Assume  < F2  . We  wish  to  show  (x)  is  convex  in 

0 £x  <_  F^(l)  . First  assume  F^  and  F2  are  absolutely  continuous. 


Then  we  need  only  show  Tllj,  (x)1  - f 

2 L 1 J 0 

nondecreasing  in  0 x <_  F^(l)  . Now 


[VI 


F2(u)du  is 


^[V0] 

— r 

dx  J 


F2  (u)du ! 


1 - V (x> 


f2|F21  V(X) 


5\« 


Let  x = iC^t)  so  that  = — p -~yt- — =r  and 

^ ' £iV«>] 


EiV(t>] 


dt  , 1L  1 vt 
dx  1 - t 


I t-Rp  (x)  1 “ 


Hence 


dH  (x) 
F1 


at  2 d F%(,)]  . 

T~  “ — - u — r~s  - and  — J F„ 

dx  1 - R_  (x)  dx  * 2 


But  F^  < F2  implies 


f,  (x) 


s’!  Fi  • ~rr- i 


f2p2  r1(») 


(u)du 


fi[rl\(>) 

£2 
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is  nondecreasing  in  0 <_  x f.F^Cl)  . Since  f”1!^  (x)  ■ t is  nondecreasing 

in  0 x <^F^(1)  , a change  of  variable  completes  the  argument. 

Since  continuous  distributions  can  be  approximated  arbitrarily  closely 
by  absolutely  continuous  distributions,  the  proof  of  part  (a)  is  complete. 

To  prove  (b)  we  will  need  the  following  fundamental  lemma. 


2.3:  Fundamental  Lemma 


If  R(0)  - 0 , 

1 x 

— J N(u)du  , then 


(2.1) 


is  nondecreasing  in  x ^ 0 and  0 <_  N(x)  <_ 


H N(u)dR(u) 
x 

J N(u)du 
0 


is  nondecreasing  in  x ^ 0 . [Note  that  if  N(x)  is  non increasing,  the 
assumption  on  N(x)  is  automatically  satisfied]. 


Proof : 


R can  be  approximated  arbitrarily  closely  from  below  by  positive 
linear  combinations  of  simple  functions  of  the  form 


0 x < x. 


X X Xq  . 


Hence  we  need  only  verify  the  lemma  for  simple  functions.  The  general 


result  follows  from  the  Lebesgue  monotone  convergence  theorem.  For  a 
simple  function,  R 
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N (u)du 


A 

‘oN(V  + / 


N (u)du 


N(u)du 


X < X, 


x lxo 


Hence,  for  x ^ x^ 


A 

j N(u)dR(u) 

0 

x 

N(u)du 


:qN (Xq)  " J N (u)du 

0 . 

x 

J“  N(u)du 


By  assumption,  N(xQ)  / N(u)du  so  that  the  lemma  follows. 


Theorem  2.1,  Part  (b) : 

Let  R(x)  - F21F1(x)  . By  assumption,  ~~  is  nondecreasing  in 
0 <_  x <.  F^Q.)  . Let  N(u)  - F^Cu)  , x - F^1(t)  and  substitute  in  (2.1) 
to  obtain 
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Fr1^) 


f V 

0 

f'V) 

/ fl 


(u)dR(u) 


(u)du 


f:1^) 


- . 

-i  1 r - 

Let  v - F F1 (u)  = R(u)  so  that  the  numerator  becomes  J F1 (u)dR(u) 


rz\M 


/ F2 (v)dv  . It  follows  from  the  Fundamental  Lemma  that 


f"1^) 


j ?2 (v)dv 


Hp1  (t)  J F2 

r 2 _ 0 

H^(t)  F’^t) 

1 / »x 


(u)du 


is  nondecreasing  in  0 t <_  1 or  is  nondecreasing  in 


0 < * < Fx  (1)  , i.e. , 


\ 5 \ 


Hp  < Hp  DOES  NOT  IMPLY  F1  < F2 


1 2 


Let  G(x)  - 1 - e so  that  HG(x)  - x for  0 < x < 1 . It  is 
easy  to  find  examples  such  that  Hp  < but  F G ; i.e.,  F is  not 
I FRA.  Note  that  for  0 < t^  < 1 , c • -log  (1  - t^)  - t^  > 0 . 

Hence 


is  not  I FRA  since 


0 <_  x < ^ 


x > t 


1 


is  decreasing  for  x >_  t^  . 
But 


Hp1 (t)  - 


0 _<  t < 

1 t < i , 


Hp1(t) 

is  anti-starshaped;  i.e.,  — is  nonincreasing  in  0 < t < 1 so  that 

«F  t HG  * 

The  significance  of  this  example  is  that  an  anti-starshaped  total 
time  on  test  plot  is  not  necessarily  evidence  that  F is  I FRA. 


If  J xdF(x)  < « , then 
0 


J xdHp(x)  - 2^x11  - F(x)]dF(x)  . 
0 0 


Proof : 


00  1 

Since  J xdHp(x)  (u)du  , we  Integrate  the  latter  by  parts  to 


r>...  - ^ 


obtain 
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J.  X X 

J H^ujdu  = j F_1(u)du  - /t(l-  t)dF_1  (t) 


Integrate  by  parts  again  to  obtain 


t(l  - t)dF 


'1(t)  ■/p'1 


(u) (1  - 2u)du 


so  that 


x X co 

y* H^OOdt  = 2 y*(l  - u)F_1(u)du  = 2 /■  [1  - F(x)]dF(x) 


by  a change  of  variable. 


Examples: 


-fix')" 

For  F(x)  * 1 - e with  mean,  say  0 , 


!/xdv 


x)  = 1/2J 


For  the  gamma  distribution 


k k-1  -Xu 


C 1 u e 
(k  - 1 


k « 1,2,  . 


with  mean  0 * — , 
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i 


fc*  ! 


• /*■»«  - 2 ftk)  t*  ■ 


i 00 

The  numerical  relationship  between  - f xdR^x)  and  the  shape 

parameter  for  Weibull  and  gamma  distributions  is  shown  in  the  following 
table. 
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